Astronomy <fc Astrophysics manuscript no. h4637 


February 2, 2008 


(DOI: will be inserted by hand later) 





BeppoShX observation of PSR B1937+21 

L. Nicastro 1 , G. Cusumano 1 , O. Lohmer 2 , M. Kramer 3 , L. Kuiper 4 , W. Hermsen 4 , T. Mineo 1 , and W. 

Becker 5 

1 IASF-CNR, Via U. La Malfa 153, 90146 Palermo, Italy 

2 Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 

3 University of Manchester, Jodrell Bank, Macclesfield SK11 9DL, United Kingdom 

4 SRON National Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands 

5 Max-Planck-Institut fur extraterrestrische Physik, Giessenbachstrafie, 85740 Garching, Germany 

Received 27 June 2003 / Accepted 10 October 2003 

Abstract. We present the results of a BeppoSAX observation of the fastest rotating pulsar known: PSR B1937+21. 
The ~ 200 ks observation (78.5 ks MECS/34 ks LECS on-source time) allowed us to investigate with high statistical 
significance both the spectral properties and the pulse profile shape. The pulse profile is clearly double peaked at 
energies > 4 keV. Peak widths are compatible with the instrumental time resolution and the second pulse lags the 
main pulse 0.52 in phase, like is the case in the radio. In the 1.3-4 keV band we detect a ~ 45% DC component; 
conversely the 4-10 keV pulsed fraction is consistent with 100%. The on-pulse spectrum is fitted with an absorbed 
power-law of spectral index ~ 1.2, harder than that of the total flux which is ~ 1.9. The total unabsorbed (2-10 
keV) flux is F 2 _io = 4.1 x 10~ 13 erg cm~ 2 s~\ implying a luminosity of L x = 5.0 x 10 31 6 (d/3.6 kpc) 2 erg s _1 
and a X-ray efficiency of r\ = 4.5 x 10 -5 O, where O is the solid angle spanned by the emission beam. These results 
are in agreement with those obtained by ASCA and a more recent Rossi-XTE observation. The hydrogen column 
density JVh ~ 2 x 10 22 cm -2 is ~ 10 times higher than expected from the radio dispersion measure and average 
Galactic density of e _ . Though it is compatible (within 2a) with the Galactic (Hi derived) value of ~ 1 x 10 22 
cm -2 , inspection of dust extinction maps reveal that the pulsar falls in a highly absorbed region. In addition, 
1.4 GHz radio map shows that the nearby (likely unrelated) Hn source 4C21.53W is part of a circular emission 
region ~ 4' across. 
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1. Introduction 

The various X-ray satellites, from ROSAT to CHANDRA, 
observed about a dozen millisecond pulsars (MSPs) (see 
e.g. Becker 2001, Becker & Aschenbach 2002). Such obser- 
vations demonstrate that the X-ray emission from MSPs 
is mainly not of thermal origin and those with the hardest 
spectra appear to be objects with strong magnetic fields 
Bl at the light-cylinder radius Rl (Rl = cP/2ir where P 
is the spin period of the MSP; see Saito et al. 1997, Kuiper 
et al. 1998, Kuiper et al. 2000, Becker 2001). However re- 
sults, in the soft "ROSAT energy band" 0.1-2.4 keV from 
a Chandra observation of 47 Tuc, suggest that significant 
thermal emission (at energies kT ~ 0.2 keV) can be pro- 
duced from MSPs (see Grindlay et al. 2002 for details). 
This appears not to be the case for PSR B1821— 24 (see 
e.g. Becker et al. 2003). Correlation between spin-down 
energy loss and X-ray luminosity was investigated by sev- 
eral authors (e.g. Verbunt et al. 1996, Becker & Triimper 
1997, Takahashi et al. 2001, Possenti et al. 2002), sug- 
gesting that the claimed law Ljc(2-10 keV) oc E 1 , with 



7 in the range 1.0-1.5, is valid for MSPs and "ordinary" 
pulsars in the same way. 

Unfortunately, good spectral and temporal informa- 
tions exist only for about half of the targeted sources. 
In fact (with the exception of the data collected by Rossi- 
XTE and, partially, by ASCA and BeppoSAX) MSP X-ray 
observations are usually affected not only by low statistics 
but sometimes also by insufficient time accuracy to per- 
form detailed periodicity and timing analyses. 

PSR B1937+21, with a period P ~ 1.56 ms, is the first 
and still the fastest rotating MSP known. Its radio pulse 
profile is double peaked with phase separation of about 
0.52. In spite of its low surface magnetic field strength of 
B s = 3.2 x\Q w (PP) 1 / 2 = 4.1xl0 8 G,its magnetic field at 
the light-cylinder is the highest of all known pulsars: Bl — 
B s (R s /Rl) 3 = 2.97 x 108PV2P-5/2 ~ 1 x 10 6 G (it is 
Rl ~ 7.4 x 10 6 and the neutron star radius Rs = 10 6 cm), 
very similar to that of the Crab pulsar. With a spin-down 
luminosity of E = Att 2 IP/P 3 ~lx 10 36 erg s" 1 (I = 10 45 
g cm 2 , momentum of inertia), the derived spin-down flux 
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density S/(47r(d/3.6 kpc) 2 ) ~ 7 x 10~ 10 erg cm" 2 s" 1 is 
similar to that of PSR B1821-24 but it is 10 times larger 
than that of PSR J0218+4232, i.e. the other two MSPs 
detected in X-rays with good statistics which have high 
luminosities and show non-thermal emission only. 

Analysis of a sensitive observation performed by Rossi- 
XTE has revealed that the X-ray peaks (see below) 
are almost perfectly aligned with the radio giant pulses 
(Cusumano et al. 2003). Such giant pulses are short 
bursts of emission with a flux density exceeding that 
of the whole profile by a factor of 10 to 100 or more. 
Having been discovered first for the Crab pulsar (Staelin 
& Reifenstein 1968), they have also been detected for 
PSR B1937+21 (Wolszczan et al. 1984, Backer 1995), and 
more recently for PSR B1821— 24 (Romani & Johnston 
2001), PSR B0540-69 (Johnston & Romani 2003) and 
PSR B1112+50 (Ershov & Kuzmin 2003). The observa- 
tions suggest that the radio giant pulses are related to the 
high energy emission of these pulsars, possibly originating 
in outer gaps (Romani & Johnston 2001). Interestingly, 
Kinkhabwala & Thorsett (2000) report that the giant 
pulses of PSR B1937+21 appear in windows of ~ 55-^70 /us 
located after the main and interpulse peaks, with a typical 
duration of less than 10 /zs. Thereby, the location of the X- 
ray peak coincides with the giant pulses and not with the 
normal radio profile, strongly supporting the idea about a 
common origin of giant pulses and high energy emission. 

In this paper we present the results of the tempo- 
ral and spectral analysis of a BeppoSAX observation 
of PSR B1937+21. Compared to those performed with 
ASCA (Takahashi et al. 2001) our observation exploits 
the better instrumental sensitivity and a longer exposure 
time. This allows the detection of the pulse profile inter- 
pulse and source photons down to energies of 0.5 keV and 
enables us to put constraints on the Na toward the pul- 
sar. Our results are compared with those obtained with 
ASCA and, more recently, with Rossi-XTE (Cusumano et 
al. 2003). 



Table 1. JPL DE200 ephemeris of PSR B1937+21 and 
derived parameters. 



Parameter 


Value 


Right Ascension (J2000) 


19 h 39 m 38! 5600084(8) 


Declination (J2000) 


21° 34' 59'.'13548(14) 


Frequency (/, Hz) 


641.928246349481(9) 


Freq. derivative (/, Hz s _1 ) 


-4.330999(2) x 10" 14 


Freq. 2 nd deriv. (/, Hz s~ 2 ) 


1.558(8) x 10~ 26 


DM a (cm" 3 pc) 


71.03998(6) 


Epoch of ephemeris (MJD) 


52328.0 


Range of validity (MJD) 


45986 - 52764 


RA proper motion (mas yr -1 ) 


-0.128(8) 


Dec proper motion (mas yr _1 ) 


-0.486(12) 



Gal latitude (I) 


57.5215 


Gal longitude (b) 


-0.2697 


Parallax (n, mas) 


< 0.28 


Distance 15 (d - from DM, kpc) 


3.6 


Surface magnetic field (Bs, G) 


4.1 x 10 8 


Light-cylinder B (B L , G) 


1.0 x 10 6 


Spin-down energy loss (E, erg s _1 ) 
E/4nd 2 (erg cm~ 2 s _1 ) 


1.1 x 10 36 


7.1 x 10~ 10 



Note: obtained from observations with the Effelsberg 
100-m telescope and Arecibo archival data. Figures in 
parentheses represent uncertainties (2a) in the last digits 
quoted. 

a Dispersion measure (DM) varies with time. 
b From NE2001 model (Cordes & Lazio 2002). 



cal background is ~ 12% higher for both LECS and 
MECS; background photons were 50 and 123, respec- 
tively. 

— a Maximum Likelihood (ML) approach to extract the 
number of counts assigned to the pulsar taking into 
account the presence of other sources and the back- 
ground simultaneously (see e.g. Kuiper et al. 1998). 



2. Observation and spatial analysis 

BeppoSAX observed PSR B1937+21 on May 1 st 2001. 
Pointing coordinates were those from the radio timing. 
The ~ 200 ks observation time resulted in a total effec- 
tive on-source time of ~ 78.5 ks and ~ 34 ks with the 
MECS (1.3-10 keV; Boella et al. 1997) and LECS (0.1- 
10 keV; Parmar et al. 1997), respectively. To extract the 
source photons from LECS and MECS we followed two 
alternative ways: 

— extraction of the photons from the circular region and 
in the energy range which optimize the signal to noise 
ratio: they were 3', 0.5-8 keV for the LECS and 2', 
1.3-10 keV for the MECS, resulting in 106 and 385 col- 
lected photons, respectively. The local background (in 
a circular region in the field of view 10' away from the 
pulsar) was compared to that obtained using archival 
data of blank sky observations. We found that the lo- 



3. Temporal analysis 

The arrival times of the extracted photons were converted 
to the Solar System Barycentric (SSB) frame using the 
BARYCONV 1 code and then searched for periodicities by 
folding the MECS data with trial frequencies around the 
radio one 2 . The timing ephemeris listed in Table ^was ob- 
tained from observations with the Effelsberg 100-m radio 
telescope and archival Arecibo data. Details are reported 
in Cusumano et al. (2003). We found a clear signal, but 
the corresponding frequency fx — 641.92825315 Hz devi- 
ated from the radio value by (/r — fx) — —7.6 x 10 -6 
Hz. This difference is about ten times the statistical error 
(~ 8 x 10 -7 Hz). It could be ascribed to incorrect pho- 
ton time markers, likely due to long term degradation of 

see www.asdc.asi.it/bepposax/software/saxdas/ 
2 Folding of the lower statistics LECS data did not reveal 
any modulation at > 2a level. 
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Fig. 1. Radio pulse profile and the 1.3-10 keV BeppoSAX 
X-ray profile phase aligned with the radio giant pulses 
(see text) marked by the two vertical dotted lines (see 
Cusumano et al. 2003). Pulse 2 lags pulse 1 in phase by 
0.52, like for the radio. The small vertical bars under the 
X-ray profile mark the phase bins used for on-pulse spec- 
tral analysis (see text). 



Fig. 2. The X-ray pulse profile in the full MECS X-ray 
band 1.3-10 keV and two sub-bands. The shaded areas 
show the estimated DC level (±lcr); the dashed lines in- 
dicate the measured background level. In the 1.3-4 keV 
energy band the secondary peak is detected at < 2a level 
above DC level; in the 4-10 keV band, instead, the detec- 
tion level is > 4cr (above background or DC level). 



the on-board clock. This effect is relevant for long dura- 
tion observations only (see Nicastro et al. 2002 for more 
details). 

Figure ^ shows the folded light curve of the MECS 
data together with the radio one obtained with the 100- 
m Effelsberg radio-telescope. The X-ray profile is dou- 
ble peaked: the primary peak (PI) is detected at 10<r 
level, the secondary peak (P2) at > 4cr. We performed a 
Gaussian fit on the two peaks. The derived phase widths 
are opi = 0.030 and <7p2 — 0.024, respectively. This values 
translate in ~ 100 /is FWHM, fully compatible with the 
instrument time resolution limit. The phase separation of 
the two peaks (primary to secondary) is 0.52 ± 0.02 (90% 
confidence), where the error is derived from the statisti- 
cal uncertainties in the peak position gaussian fits: ±0.007 
and ±0.010 for PI and P2, respectively, Figure El shows 
the X-ray pulse profile in 3 energy ranges. Non-detection 
of the interpulse in the ASCA data is explained by their 
lower statistical significance (227 total photons in 1.7-6.5 
keV). 

To quantify the pulsar DC component, we applied the 
bootstrap method proposed by Swanepoel et al. (1996), 
which allows us to estimate the off-pulse level (sky back- 
ground plus unpulscd source component) in the pulse pro- 
file. We find that over the total energy range 1.3-10 keV 



the source has a pulsed fraction 3 of / = 85 ± 5% reduc- 
ing to 54 ± 7% at energies 1.3-4 keV. In the 4-10 keV 
band the detected photons are ~ 100% pulsed. Changing 
the folding initial phase causes the pulse profile to change 
slightly but the detection level of the two peaks and the 
pulsed fraction remain well within the statistical error. 
The ASCA pulsed fraction of 59% is compatible with our 
results. 

There is an indication for an increase of the strength 
of P2 compared to PI with energy. Using the bins marked 
in Fig.^as on-pulse phases, P1/P2 = 5.13 ± 1.69 for the 
band 1.3-4 keV and P1/P2 = 2.99 ± 0.84 for the band 
4-10 keV. 

A joint analysis of radio data (Jodrell Bank) and X- 
ray data (ASCA) led Takahashi et al. (2001) to claim an 
alignment of the main X-ray peak with the radio inter- 
pulse. The BeppoSAX on-board clock absolute accuracy 
does not allow us to align our profile to the radio one. 
However we noticed earlier (Nicastro et al. 2002) that the 
radio and X-ray peaks have both a phase separation of 
~ 0.52 (primary to secondary), then leading us to the con- 
clusion that their result is likely incorrect. Results from a 

3 It is / = the percentage of total counts above DC level 
with respect to the net total counts in the full profile above 
the background level (determined in the spatial analysis). 
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Table 2. Power-law fit parameters. 



Spectrum 




a 


N H 
(f0 22 cm- 2 ) 


^2-10 


Total 


1 5+ a2 


1 94+ - 13 

■■■•"-0.11 


9 i R+0.90 
z - lo -0.65 


4.1 ±0.6 


Pulsed c 


0.33 


-, O1 +0.15 

l - i± —o.i3 


2 1+ 19 
z - 1 -i.i 


3.0 ±0.5 



Note: all quoted uncertainties are 2a (95%) confidence. 

Xn — 0.96 for both the Total and Pulsed spectra. Photons 

distribution law N p h{E) = K x E~ a . 
a Normalization at 1 keV (10~ 4 ph s _1 cm -2 ). 
b Unabsorbed flux in 2-10 keV (10" 13 erg cm -2 s _1 ). 
c MECS 1.3-10 keV (see marked bins in Pig. [Q. 



Rossi-XTE observation has confirmed this doubt showing 
that the X-ray peaks are almost perfectly aligned with the 
radio giant pulses (see Cusumano et al. 2003). 

4. Spectral analysis 

In spite of the lower exposure time and sensitivity of LECS 
we were able to detect a significant "unpulsed" signal in 
the 0.5-8 keV range for this instrument. Inclusion of these 
softer photons in the spectral fit allowed us to better con- 
strain the 2Vh toward the pulsar. In grouping the LECS 
and MECS data we applied the condition to have at least 
20 and 30 counts per bin, respectively. The best fitting 
model (in the detection range 0.5-10 keV) was an ab- 
sorbed power-law. A Black Body model did not fit the 
total flux data satisfactorily (kT ~ 1.2 keV and Xn — 1-5, 
13 d.o.f.). We then tried to fit the data with a broken 
power-law and a power-law + Black Body. For the latter 
case we also tried fixing the value of kT to various values 
in the range 0.1 ~ 0.5 keV as higher values are unlikely for 
an MSP (Zavlin et al. (2002) found kT ~ 0.2 keV for PSR 
J0437-4715 and Grindlay et al. (2002) report kT ~ 0.22 
keV for the MSPs in 47 Tuc). The F-test gives chance im- 
provement probabilities of 90% for the broken power-law 
model and values around 50% for the power-law ± Black 
Body, then we can reject these models. In the energy bands 
0.1-2.4 and 2-10 keV we can estimate absorbed flux upper 
limits of ~ 1 x 10~ 14 and ~5x 10~ 14 erg cm -2 s _1 for 
a Black Body component with kT — 0.5 keV. On-pulsc 
analysis, subtracting the off-pulse counts, was also per- 
formed; phase bins used to extract the pulsed photons are 
marked in Fig.Q] Fit of the 1.3-10 keV MECS data alone 
gave Nh = (1.2 ± 0.9) x 10 22 cm" 2 and a = 1.7 ± 0.4 
(X 2 „ = 1-06). 

In Table [5] the best fit results (for the absorbed 
power-law model) are reported; in Fig. [21 the total 
(LECS±MECS) and pulsed spectra are shown together 
with the unabsorbed power-law and the fit residuals for 
both the total and the pulsed spectra. We also fitted the 
photons collected in the off-pulse region (DC level) which 
gave, though at lower statistics, a consistent Nh and a 
photon index of ~ 3.3 ± 0.5. We note that the spectral 
index of 1.2 of the pulsed emission is harder than that 




l 3 7 10 

Energy (keV) 



Fig. 3. The 1.3-10 keV MECS (open diamonds) plus 0.5- 
8 keV LECS (open squares) total spectrum together with 
the MECS pulsed spectrum (filled circles). The straight 
line is the unabsorbed power-law with a = 1.94. The 
strong absorption below ~ 3 keV is evident. Fit residu- 
als are also shown. 



obtained from the total emission. This is a typical value 
for magnetospheric emission and is fully consistent with 
the value obtained by Cusumano et al. (2003) using Rossi- 
XTE data. The total unabsorbed and absorbed fluxes of 
PSR B1937±21 in the 2-10 keV range are 4.1 x 10" 13 
and 3.5 x 10 -13 erg cm -2 s _1 , respectively. The derived 
(unabsorbed) luminosity is Lx = 5.0 x 10 31 8 (d/3.6 
kpc) 2 erg s -1 while the X-ray conversion efficiency is 
j] = Lx/E = 4.5 x 10~ 5 0, where 6 is the solid angle 
spanned by the emission beam. This result is confirmed 
by the Rossi-XTE 2.5-17 keV spectrum (Cusumano et al. 
2003). We also note that the fit parameters are consis- 
tent with the non-detection (Verbunt et al. 1996) in the 
ROSAT HRI energy range (0.1-2.4 keV). 

5. Discussion 

5.1. Neutral hydrogen column density 

Galactic Nh can be estimated from radio/IR measure- 
ments in three different ways: 1. from the DM by assum- 
ing 10 H atoms for each e _ ; 2. from Hi measurements, 
assuming there is no hydrogen at molecular state; 3. from 
dust extinction using the relation Nh = 1.79 x 10 21 Av of 
Predehl & Schmitt (1995). When there are no direct pul- 
sar measurements or three dimensional distribution mod- 
els, points 2. and 3. make use of models of the integral 
content of H along a particular line of sight. They then 
only give upper limits for the atomic and total hydrogen 
content, respectively. 
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For PSR B1821-24 placed in the globular cluster M28, 
at a distance of ~ 5.5 kpc and galactic latitude of — 5?6, 
one finds TVjj — 1.6 x 10 21 cm^ 2 which is consistent with 
the estimates from Galactic H I, dust and pulsar dispersion 
measure (Becker et al. 2003). Similarly consistent results 
are obtained for PSR J0218+4232. However, because of 
variable e~ /iVn ratios in the various Galactic directions 
and the presence of more o less dense gas clouds along the 
line of sight, discrepancies would not be surprising. 

The column density we derive from the power-law fit 
is 10 times greater than the radio DM derived one and is 
only marginally consistent (within 2a) with the H I derived 
Galactic value 1.2 x 10 22 cm~ 2 (Dickey & Lockman 1990). 
In addition, it is apparently totally inconsistent with the 
7.4 x 10 22 cm" 2 obtained from the Schlegel, Finkbeiner & 
Davis (1998; SFD) dust extinction maps (Schlegel et al. 
(1998). 

Being the line of sight to PSR B1937+21 tangent to the 
Galactic spiral arms (/ ~ 57.51, b ~ —0.29, see e.g. Fig. 
3 of Watson et al. 2003), one can invoke a particularly 
low e^/A^H ratio in that direction due to the lack of ioniz- 
ing sources. H I maps report the integral Galactic column 
density but do not account for hydrogen in molecular form 
and then, in particular for objects in the Galactic plane, 
they must generally be regarded as lower limits. A more 
reliable check of the X-ray determined A^h can come from 
the investigation of dust extinction maps. Using the afore 
mentioned Predehl & Schmitt (1995) relation, Verbunt at 
al. (1996) set an upper limit of 2.2 x 10 21 cm -2 for the 
column density (i.e. in agreement with the DM value). 
From the more recent SFD map, as shown in Fig. we 
note that PSR B1937+21 is almost at the center of an en- 
hanced absorption region. The map has a pixel resolution 
of 2(37 x 2(37 and spatial resolution of 6.'1 (FWHM). The 
enhanced region is then consistent with a point source at 
coordinates RA (J2000) ~ 19 h 39 m 5 Dec (J2000) ~ 21° 
37.'5. Given the mentioned spatial resolution, it is com- 
patible with the position of PSR B1937+21. However at 
this position there is a strong radio and IR source: IRAS 
19375+2130 (= 4C21.53W = G57-0.27 = G57.55-0.27). 
This source was extensively investigated since the discov- 
ery of PSR B1937+21 (see e.g. Becker & Helfand 1983). 
It is listed as an ultra-compact Hn region and selected 
by Wood & Churchwell (1989) as a candidate site of OB 
star formation. To estimate the visual Galactic extinction 
toward the pulsar, we used the map value of neighboring 
regions: it is E(B— V) ~ 6 and a value of Ay ~ 19 is 
derived, giving ATh = 3.3 x 10 22 cm -2 . This is some 50% 
higher than the X-ray derived N H for PSR B1937+21. The 
observed high dust/gas ratio can then be justified with- 
out requiring a very large distance to the pulsar. However 
it has to be larger than the NE2001 model derived one 
(Cordes & Lazio 2002) as 3.6 kpc is a lower limit de- 
rived from timing (see below). Considering that the SFD 
map predictions are not very reliable at galactic latitudes 
|6| < 5° (Schlegel et al. 1998), the consistency with the 
X-ray derived value is remarkable. 



5.2. The pulsar distance 

Timing observations of the pulsar provide proper motion 
measurements and an upper limit on the timing parallax of 
7r < 0.28 mas (see Table 0. This information can be used 
to infer that the DM distance of 3.6 kpc is, in fact, a lower 
limit. An upper limit can be derived from the transversal 
velocity estimate given by the interstellar scintillation ve- 
locity viss ~ 50 ± 10 km s" 1 found by Cordes (1990). 
Anyway, given the timing derived pulsar proper motion, 
the resulting limit of d = u/(4.74 fi) = u/(4.74 x 0.5) = 21 
kpc is not very constraining. While it has been shown that 
the scintillation velocity is indeed a good estimator for the 
true transversal velocity (Nicastro et al. 2001), for low ve- 
locities such as measured for MSPs (i.e. ~ 30 km s _1 ), 
systematic effects can result in large errors. A campaign 
of scintillation study could lead to a more reliable viss es- 
timate, but it is unlikely that a value lower than 30 km 
s _1 can be measured and therefore it can only be useful to 
set the Galactic upper limit distance of < 10 kpc. On the 
other hand, regular multi-frequency observations can help 
to remove the effects of "interstellar weather" , which de- 
grades the timing precision of PSR B1937+21 significantly 
(e.g. Backer & Wong 1996). If these effects can be mod- 
eled, higher timing precision is obtainable in the future, 
possibly helping in setting tighter limits on the parallax 
measurements. 

Heiles et al. (1983) studied the Hi absorption spectra 
toward PSR B1937+21, 4C21.53W and another close-by 
source. By setting the Galactic tangent point to 5.4 kpc, a 
comparison of position and intensity of the lines for these 
three objects allowed them to conclude that 4C21.53W is 
at ~ 10.7 kpc and the pulsar at ~ 5 kpc from us. Frail & 
Weisberg (1990) revised this result using a new model for 
the rotation of the Galaxy. They adopt the tangent point 
distance as the pulsar distance lower limit: d > 4.6 ± 1.9 
kpc. Watson et al. (2003) have recently confirmed the po- 
sition of the H II region in the Perseus arm. Assuming a 
tangent point at a distance of 4.6 kpc, they derive a dis- 
tance of 8.6 kpc by studying H2CO absorption and HllOa 
emission lines (HllOa is the 4874.1570 MHz radio recom- 
bination line). However the data are inconclusive on the 
position of the intervening molecular cloud, giving either 
3.6 or 5.5 kpc as possible distances, for the two sides of 
the tangent point, respectively. With the pulsar at a dis- 
tance > 3.6 kpc, it is possible that the cloud contributes 
significantly to the absorption of the low energy X-rays. 

We also investigated the 1.4 GHz NVSS data (Condon 
et al. 1998). Figure shows the 15' x 15' region around 
the pulsar position. PSR B1937+21 lies on the edge of a 
diffuse circular region (flux is ~ 50 4- 60 mJy/beam) and 
is moving away from it as indicated by the proper motion 
vector shown in Fig. Incidentally, the vector orientation 
along the Galactic plane can explain why PSR B1937+21 
has the lowest z-height of all MSPs. The strong source 
to the North is the aforementioned 4C21.53W and has a 
flux density of ~ 800 mJy, compared to the ~ 20 mJy 
of PSR B1937+21. The non-detection of the diffuse re- 
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Fig. 4. 8° x 8° SFD (Schlegel et al. 1998) dust extinc- 
tion map centered on PSR B1937+21. The Galactic plane 
is clearly visible. The circle (24' diameter) is centered on 
IRAS 19375+2130. The extinction toward this source is 
about three times higher than toward the surrounding re- 
gions ~ 20' away from it. The pulsar position is marked 
with a cross in the zoomed inset. 



Fig. 5. 0?25 x 0?25 NVSS 1.4 GHz image around 
PSR B1937+21. The pulsar is marked with a cross with 
the proper motion direction toward South- West shown. 
The diffuse region is ~ 4' across with ~ 60 mJy/beam 
at the center. The strong source is 4C21.53W (~ 360 
mJy/beam). Contours are log spaced. 



gion in the VLA maps by Becker & Helfand (1983) can be 
explained by their use of a different array configuration. 
They used the B-array rather than the D-configuration as 
used for the NVSS, resulting in a much degraded sensitiv- 
ity for diffuse emission. Though the circular shape of the 
region and the pulsar velocity vector could be suggestive 
of an association, the curved (externally "compressed") 
shape of the Hn region indicates a relation with the lat- 
ter (see Fig. 1 of Becker & Helfand 1983). A possible way 
to estimate the chance coincidence of PSR B1937+21 with 
the H ii "complex" is to select the sources of Watson et 
al. (2003) away from the inner arms, i.e. 70° < I < 50° 
and |6| < 1°, and assume that they all have an associated 
"region" of 4' diameter (we recall that those Hn regions 
were selected from the IRAS point source catalogue, i.e. 
they have FWHM < 2' at 100 fxm). The ratio of the total 
area of the selected 20 sources to 20 x 2 square degrees 
is 9 x 10~ 4 . We note also that in the same area there are 
15 known pulsars, two of which are MSPs 4 . The coinci- 
dence at this galactic longitude of a peculiar pulsar with 
such a small region (at a distance of 8.6 kpc the region 
would be ~ 20 pc across, ~ 3.5 pc at 3.6 kpc) remains 
then remarkable. 

5.3. Spectral index and magnetospheric emission 

Table lists the main emission properties of the three 
MSPs which show non-thermal emission only. In order 
to account for the contribution to the DC flux by the 

4 http:/ /www. atnf.csiro.au/research/pulsar/psrcat/ 



field source listed by Becker et al. (2003), but unresolved 
by the BeppoSAX MECS, we re- analysed the BeppoSAX 
observation of B1821— 24. The pulsed fraction was com- 
patible with 100%. This is consistent with the result by 
Danner et al. (1997) using ROSAT HRI data in the 0.1-2 
keV band and explains the similar indices of the on-pulse 
BeppoSAX spectrum and the total CHANDRA spectrum, 
for which no pulse phase selection could be performed. As 
can be seen from Table the "total" flux photon index of 
PSR B1937+21 is significantly softer than that measured 
for PSR J0218+4232 and PSR B1821-24. Also, we found 
that the spectrum of the pointlike DC component of PSR 
B1937+21 is softer than its pulsed spectrum. It is interest- 
ing to note that high-resolution CHANDRA observations 
of PSR J0218+4232 (Kuiper et al. 2002) confirmed the 
presence of a DC component of this source, earlier noted 
in Kuiper et al. (1998) and Mineo et al. (2000). The spec- 
trum of the DC component of PSR J0218+4232 appeared 
to be softer than its pulsed component, like we found here 
for PSR B1937+21. Kuiper et al. (2000) concluded that 
we do not expect thermal emission from the surface of an 
old, recycled MSP due to cooling, and that the soft DC 
component is most likely due to reheating of the polar 
cap area by backfl owing accelerated particles. This has to 
be investigated further for PSR B1937+21. A high spa- 
tial resolution observation may reveal field sources and/or 
diffuse emission also in the PSR B1937+21 field. Such in- 
formation would allow us to better constrain both its to- 
tal and pulsed emission properties as compared to those 
of PSR B1821-24 and PSR J0218+4232 for which higher 
resolution observations already exist. In spite of the large 
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Table 3. Energetics of PSR B1937+21 compared to 
B1821-24 and J0218+4232. 



Parameter 


"D 1 t 11 a 

-tSlyo / +ZL 




JU218 


Period (P, ms) 


1.558 


3.054 


2.323 


P (xl(T 19 s Hz) 
E (xlO 36 erg s" 1 ) 


1.0 


16 


0.77 


1.1 


2.2 


0.25 


B L (xlO 6 G) 


1.0 


0.77 


0.32 


Distance (d, kpc) 
E/{4ird 2 ) d (xKT 10 ) 


3.6 (4.6) 


5.5 


5.7 (4.0) 


7.1 (4.3) 


6.1 


0.64 (1.3) 


a total 


1.9412;!? 




n 04+0.22 

"•"-0.22 


a pulsed 


_. _ _ _i_n i k 
1 - zl -0.13 


i -- LU -().2<) 


0.6ll°;i 


Pulsed frac. c (/, %) 


86 


98 


73 


L { X , P (xlO 32 ) 


0.37 (0.60) 


0.62 


1.3 (0.64) 


Lx,i 


6.4 (10) 


12 


16 (7.9) 


Lx,E 


5.7 


16 


0.67 


Lx,p/E (xlO" 4 ) 


0.34 (0.55) 


0.28 


4.8 (2.4) 


L x ,i/E (xlO" 4 ) 


5.8 (9.5) 


5.6 


65 (32) 



Note: luminosity Lx (erg cm 2 s 1 ) and flux Fx (erg s x ) 

in 2-10 keV. 
a In parenthesis figures for d — 4.6 kpc. 
b See Becker et al. (2003); Kawai & Saito (1999). 
c See Mineo et al. (2000). In parenthesis figures for d — 4.0 

kpc (Bassa et al. 2003). Using Rossi-XTE, Kuiper et al. 

(2003) find over the range 2 to ~ 20 keV 1.141^ for a 

pulsed. 

d Spin-down flux (erg cm~ 2 s _1 ). 

In 1.6-10 keV as derived by BeppoSAX observations. 

f Luminosity: P = pulsed (d 2 xFx,p), I = isotropic total, E = 

isotropic expected from the relation logLx = 1-45 log E — 

19.5 (Takahashi et al. 2001). 



P difference between PSR B1937+21 and B1821-24, the 
X-ray production efficiencies of these two pulsars are sim- 
ilar; and similar are their Bl and pulsed emission spec- 
tral indices. PSR J0218+4232 is a factor of 10 more effi- 
cient even though its B>l is the weakest of the three; on 
the other hand it has the hardest spectrum. We finally 
note that using 4.6 kpc as the distance to PSR B1937+21 
and as the distance to PSR J0218+4232 the 4 kpc up- 
per limit recently set by Bassa et al. (2003), the "pulsed" 
X-ray luminosities of the three pulsars are almost iden- 
tical. Identical are also the pulsed spectral indices if for 
PSR J0218+4232 we adopt the statistically more signifi- 
cant value of a = 1.14 from Kuiper et al. (2003). 

All together these figures are not inconsistent with the 
claiming of a single Lx ~ E law for MSPs and ordinary 
pulsars suggesting that a common magnetospheric emis- 
sion process could be responsible for most, if not all, of the 
emitted X-ray flux in the band 2-10 keV. It is not clear 
if pulsar in 47 Tuc could be peculiar in this respect given 
that they have been studied only in the soft energy range 
0.1-2.4 keV (Grindlay et al. 2002). 



6. Conclusions 

We detected the double peak profile of the fastest ro- 
tating pulsar known. The 1.3-10 keV pulsed fraction is 
85% which becomes 100% in the 4-10 keV band. The DC 
component is significantly detected only in the low energy 
band 1.3-4 keV (46 ± 7% unpulscd photons). The pulse 
phase resolved spectral analysis confirms this result. The 
secondary (X-ray) peak is detected at high significance 
above 3-4 keV and the ratio primary /secondary decreases 
with energy. This suggests that the secondary peak has a 
harder spectrum. It was not possible to perform an ab- 
solute phase comparison with the radio profile, but phase 
separation comparison indicates, contrarily to the ASCA 
finding by Takahashi et al (2001), that the main radio and 
X-ray peaks are aligned. This indeed was confirmed by a 
Rossi-XTE observation which shows that the X-ray peaks 
are almost perfectly aligned with the radio giant pulses 
(Cusumano et al. 2003). 

We measure a hard power-law spectral index a ~ 1.2 
for the pulsed photons. This is similar to the value found 
for PSR B1821-24 and PSR J0218+4232 (Kuiper et al. 
2003). We find a column density of iV H ~ 2.1 x 10 22 
cm~ 2 toward PSR B1937+21. It is marginally consistent 
with the integral Galactic value obtained by the Dickey & 
Lockman (1990) Hi model of 1.2 x 10 22 cm~ 2 but in agree- 
ment with the SFD map prediction, and gives a e~/Nu 
ratio of 1/100. 

From the NVSS maps we note that the pulsar is located 
on the edge of a diffuse ~ 4' circular region with the strong 
radio/IR source 4C21.53W and the pulsar located at op- 
posite edges. The pulsar also seems to move away from this 
region with a transverse speed vector along the Galactic 
plane. Emission/absorption line observations have shown 
that 4C21.53W is very likely located in the Perseus arm 
at 8-9 kpc away from us (Watson et al. 2003) whereas it 
is unlikely the pulsar to be too far away from the Galactic 
tangent point, i.e. 4.6 kpc (Frail & Weisberg 1990). These 
results and the morphology of 4C21.53W/diffuse region 
suggest they are associated, with the pulsar aligned by 
chance. 

A detailed study of the ISM toward the pulsar could 
possibly help modeling the (relatively) high timing noise 
of PSR B1937+21 and consequently further constrain the 
parallax mesurement. However, we consider ~ 5 kpc to be 
a good distance estimate. 

A high spatial resolution CHANDRA X-ray observa- 
tion can help both to verify the high A% value and to 
study the continuum X-ray emission properties of its close 
surroundings. A better resolved 1.4 GHz radio map (ob- 
tainable from C configuration VLA archival data) can help 
to study the morphology of the "4C21.53 complex" and 
to confirm the chance alignment of the diffuse emission 
region with the pulsar. 
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